Tunneling spectroscopy of Luttinger-liquid structures far from equilibrium 
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We develop a theory ol tunneling spectroscopy of interacting electrons in a non-equilibrium quan- 
tum wire coupled to reservoirs. The problem is modelled as an out-of-equilibrium Luttinger liquid 
with spatially dependent interaction. The interaction leads to the renormalization of the tunneling 
density of states, as well as to the redistribution of electrons over energies. Energy relaxation is con- 
trolled by plasmon scattering at the boundaries between regions with different interaction strength, 
and affects the distribution function of electrons in the wire as well as that of electrons emitted from 
the interacting regions into non-interacting electrodes. 

PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 
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I. INTRODUCTION 

One-dimensional (ID) interacting fermionic systems 
show remarkable physical properties and are promis- 
ing elements for future nanoelectronics. The electron- 
electron interaction manifests itself in a particularly dra- 
matic way in ID systems, inducing a strongly correlated 
electronic state - Luttinger liquid (LL)ii2jiLi. A paradig- 
matic experimental realization of quantum wires are car- 
bon nanotubes^; for a recent review see Ref. H. Further 
realizations encompass semiconductor, metallic and poly- 
mer nanowires, as well as quantum Hall edges. 

There is currently a growing interest in non- 
equilibrium phenomena on nanoscalcs. A tunneling spec- 
troscopy (TS) technique for non-equilibrium nanostruc- 
tures was developed in Ref. 7. Employing a supercon- 
ducting tunneling electrode allows one to explore not only 
the tunneling density of states (TDOS) but also the en- 
ergy distribution function. The energy relaxation found 
in this way provides information about inelastic scatter- 
ing in the system. In a very recent experiment^ this TS 
method was applied to a carbon nanotube under strongly 
non-equilibrium conditions. 

In this paper, we develop a theory of TS of a LL out 
of equilibrium. Specifically, we consider a LL conductor 
connected, via non-interacting leads, to reservoirs with 
different electrochemical potentials, (il — = eV and 
different temperatures Tl, Tr (where the indices L, R 
stand for left- and right-movers). It is assumed that the 
coupling to the leads is adiabatic on the scale of the Fermi 
wave length, so that no backscattering of electrons takes 
place. We model the leads as non-interacting ID wires, 
so that the electron-electron interaction is turned on at 
the vicinity of the points x = ±L/2, see Fig. [T] This 
model is quite generic to properly describe the problem 
at hand, independently of the actual geometry of the 
leads. Note also that the ID setup with strongly non- 
uniform interaction may be experimentally realized by 
using external screening gates. 



It is known that energy relaxation is absent in a uni- 
form clean LL. Within the golden-rule framework, the 
lack of energy relaxation for forward scattering processes 
results from ID kinematic constraints that do not allow 
to satisfy the energy and momentum conservation laws 
simultaneously 9 . On a more formal level, the conserva- 
tion of energies of individual particles in a spatially uni- 
form LL is protected by the integrability of the system, 
which implies an infinite number of conservation laws**. 
Inclusion of spatial dependence into the model violates 
these laws and leads to energy relaxation that takes place 
at the regions where the interaction varies in spaced. 

The fact that inhomogeneous interaction induces en- 
ergy relaxation of electrons has been pointed out for the 
first time in Ref. [13| in the context of interacting quantum 
Hall edges but a detailed analysis of this effect has been 
missing until now. On the other hand, one may expect 
this to be a dominant effect on the electron distribution 
function in experiments done on modern high-quality 
quantum wires (such as ultraclean carbon nanotubes 1 ^) , 
under non-equilibrium conditions. There is thus a clear 
need in the theory of TS in non-equilibrium LL. 

It is worth noting that we assume the absence of 
backscattering due to impurities in the wire. When 
present, such impurities strongly affect the electronic 
properties of a LL wire: they induce diffusive dynamics 
at sufficiently high temperature T and localization phe- 
nomena proliferating with lowering T (Ref. [lH ) , as well as 
inelastic processe a 16 d 7 . We also neglect the nonlinearity 
of the electron dispersion whose influence on spectral and 
kinetic prop erties of ID electrons was recently studied in 
Refs. H [nj 



II. FORMALISM 

Within the LL model, the electron field is decoupled 
in a sum of right- and left-moving terms, ip(x,t) = 
ipn(x, t)e lppX +iPl(x, t)e~ lpFX , wherepF is the Fermi mo- 
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FIG. 1: Schematic view of a LL conductor with various posi- 
tions of tunnel probes. The solid curve in the lower part of the 
figure shows a spatially dependent LL interaction parameter 
K(x). The dashed line corresponds to the limit of a sharp 
variation of K(x) at the boundaries. 



mentum. The Hamiltonian of the system reads 

H = H Q +H int , (1) 
H a = -iv / dx Oip R d x ipR - , (2) 



dxg(x)(ip R i> R + iplipL) 



(3) 



where v is the electron velocity and g{x) is the spatially 
dependent electron-electron interaction constant. 

We will proceed by following the lines of the func- 
tional bosonization approach 18 in the non-equilibrium 
(Keldysh) formulatio n 19 i 20 ■ 16 . Performing the Hubbard- 
Stratonovich transformation, one decouples the interac- 
tion term via a bosonic field cf> and gets the action 



rj=R,L 



^-^--^rV, (4) 



where du t L — d t ± vd x and the fields are defined on the 
Keldysh time contour. The information about physical 
observables is contained in Keldysh Green functions 21 
G > and G < ; see, in particular, Appendix [A] where we 
express tunneling current in terms of functions G^ and 
discuss how its measurement allows to determine G^ ex- 
perimentally. The Green functions G< can be presented 
in the form 



(5) 



G^(x,t;x',t') = jV<t>Z[<t>]e-*t>° 
(x,t;x',t'), 



x G|[ 



where we introduced the Green function in a given field 
configuration, G^ [4>], and the sum of vacuum loops, Z [</>]. 

In ID geometry the coupling between the fermionic 
and bosonic fields can be eliminated by a gauge transfor- 
mation ip v (x,t) — > tl>n(x,t)e ien ( x ' t if we require 



id v v = 



(6) 



As a result, G^ [4>] can be cast in the form 
G}[4>]{x,t-x\t') = Gloix-x'-t-t'y-^^/ 2 



x e 



?|(x,i;x',t') 



(7) 



Here 



§<(x,t;x ,t') = i6± iV (x,t) - i6 TtV (x',t') 



G< is the Green function of free fermions, 



2v smh'KT ri (r]^ rl ± iO) 



(8) 



(9) 



the coordinate £,r/l = x/v =F t labels the trajectory of 
a particle, and we use the convention that in formulas 
r\ should be understood as r\ = ±1 for right/left moving 
electrons. 

It is convenient to perform a rotation in Keldysh space, 
thus decomposing fields into classical and quantum com- 
ponents, 01,(^2 = (0+ ± </'-)/v / 2, where the indices + 
and — refer to the fields on two branches of the Keldysh 
contour. Further, we introduce vector notations by com- 
bining ipi and <f>2 in a 2-vector </>. To proceed further, we 
resolve Eq. ^ and express 9 V through <f> as 



where Gno is the Green function of free bosons, 



(10) 



yr >° ~ 1 g*o o 



(ii) 



Its retarded and advanced components are given by 
7)0 ~ LU-r]vq± iO 



(12) 



The Keldysh component of G-qO is given by G^q = {G^o 



G%o)Bv wnere B-n"' ( w ) ls determined by the temper- 
ature T v of the reservoir from which the electrons moving 
in direction r/ emerge, 



>(o) 



B^>{uj) = cothw/2T, ; . 



(13) 



Using Eqs. §8§ and (fTP|) and performing a transforma- 
tion to the coordinate space, we express the exponent 
<&^[4>](x, t, x',t') through the bosonic field <f>(y), 

^[4>](x,t,x',t') = J — 
The components of J are found as 

7^ 



jl n Av) = ^ [e\,{x-y)}e-^-eW-y)}e-^) , 
_1 (e^ y _ x )] e ~i^v + e[ v (y - x')]e- l <) ,(15) 



T- 



where 6[x) is the Heviside ^-function. The vacuum loop 
factor in Eq. ([5]) is given by 



Z[<t>\ =exp( n</> 



(16) 
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where II is the polarization operator, 



III. SHARP BOUNDARIES 
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It can be decomposed into left and right moving parts, 
II = U R + U L , with 



n 



1 



R,L 



2ir oj + =f vfH 



jra _ 



1 



2tt uj- =F VFq ' 



(17) 



where u>± = u> ± iO. Performing the averaging over the 
auxiliary field </>, we get 

G*(x,t;x',t') = G| (e, - Qe-^VW\^< ; (18) 

where the effect of the interaction is represented by the 
"Debye- Waller factor" e » with 



(IhJ 



T<{x,t;x ,t ) = -- / —dyidy 2 

x ^^1)^(^^2)^(2/2). (19) 



Here 



v = (n + a-Vi)- 1 



(20) 



is the screened electron-electron interaction potential. Its 
retarded component is given by 



v r M) = g{y) 



7T 



,(21) 



where the function is determined by the following 
differential equation 



{^ + d v u 2 {y)d v )gi{ y ,y')=5{ y -y') 



(22) 



which describes the plasmon propagation in a medium 
with spatially dependent sound velocity u(x) = v(l + 
g(x)/TTv) 1 / 2 . The Keldysh component of the interaction 
propagator is obtained as 



V* (yi, y 2 ) 

where 



2ttv< 



Y,B n WMHM, (23) 



I2(y)= / dy'e^y'v:(y,y'). 



(24) 



At equilibrium, B r (uj) = Bl(uj) = B(u>), this reduces to 



V 



K 



v r - v a 



B(lo) 



(25) 



in agreement with the fluctuation-dissipation theorem. 



So far we made no restriction on the way the interac- 
tion changes in space. Let us consider first the case when 
the interaction turns on and off sharply on the scale set 
by the temperatures, It ~ vj max-fT^, T R }. This limit 
can be modelled via a stepwise interaction as represented 
by the dashed line in Fig. [1] Equation (f2"2"| for can 
be then straightforwardly solved by using the fact that 
the velocity u is constant in each of three regions and 
employing the proper boundary conditions [continuity of 
Ql{y,y') and of u 2 (y)d y gi(y,y')\ at y = ±L/2. 

In the TS context, we are interested in the Green func- 
tions G< with coinciding spatial arguments, x = x' . As- 
suming x to be in the interacting part of the wire (and 
not too close to the boundaries) and setting t' = 0, we 
find 



/ duj 



(i^) 2 <H + (i + «H 



(0), 



2(1 + K 2 ) 



x ( 1 — cos tot ) ± i sin tot 



where 



K = v/u — (1 + g/irv) 



-1/2 



(26) 



(27) 



is the conventional dimcnsionlcss parameter characteriz- 
ing the interaction strength in a LL and 



7 = 



(K-lf 



2K 



(28) 



The integral in Eq. (|26|) and in analogous formulas below 
is logarithmically divergent at large frequencies and re- 
quire an ultraviolet regularization. Specifically, these in- 
tegrals are understood as regularized by a factor e~oj/A, 
where A is an ultraviolet cutoff. Deriving Eq. (|26|) . we 
have neglected terms of the form e mujL / u (with non-zero 
integer n) that arise due to the Fabry-Perot-type inter- 
ference of plasmon modes reflected at the boundaries. 
Keeping these terms would lead to an additional oscilla- 
tory structure in energy 2 ^ with the scale iru/L. Since we 
are interested in TS of long wires, we assume that this 
scale is much less than ma.x{T R , Tl}, so that oscillations 
are suppressed. 

Substituting Eq. (J26]) into Eq. ([18|), we finally get the 
Green functions: 



G R (t) = (27TWP 

where 



-1 l + Q 



-irjeVt/2 



(^ -I) 4 R {K 2 -lf 



4K(1 + K 2 ) 



4K(1 + K 2 ) 



(29) 
(30) 



The Green functions (|29|) can be determined experimen- 
tally from TS measurements^, see Appendix fXJ Their 
difference determines the TDOS v{e) 1 



G> (e, x, x) - G< (e, x, x) = -2mv v {e), 



(31) 
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while each of them separately (or their sum) contains also 
information about the distribution function, as discussed 
below. The results for the TDOS have been found in 
Ref. M- 

Next we consider the non-interacting parts of the wire, 
and discuss, e.g., the right moving electrons. In the re- 
gion I (see Fig. [I}, x , x ' < ~L/2, we find from Eqs. (fill)) . 
(fT5| that = 0, so that the Green functions of the 
right movers are not modified by interaction. Physically 
this is quite transparent: the right-moving electrons in 
this part of the system are just coming from the reser- 
voir and are not yet "aware" of the interaction with the 
left-movers. 

The situation is distinctly different in the region III, 
x,x' > L/2. Assuming x = x\ we find 



dio (1 - K) 2 , 
— '—{I - cos ujt) 



, Q cj 1+K 
Substituting Eq. ([32]) into Eq. (JTHJ), one gets 



where 



2K 
1 + K 2 



1 T r 



n 



-irieVt/2 



K : 



1 + K 2 



(32) 



(33) 



(34) 



Since T< in Eq. (32]) is real, the TDOS is not affected by 
the interaction, ^r(e) = v a = l/27ru, as expected. The 

modification of the functions as compared to that of 
incoming electrons, G\ , implies therefore the change in 
the distribution function n#(e) of right-movers. Indeed, 
for non-interacting particles G^ = 2'KiVQn^e) and G^ = 
— 2niv [l — nji(e)]. We thus see that the electrons ejected 
from the interacting part of the wire into the lead are 
affected by the interaction: their distribution function 
has changed. 

The left-moving electrons can be analyzed in the same 
way; the corresponding results are obtained by replacing 
R «-> L, V <-> -V in Eqs. (H, (33]). Clearly, the role of 
the regions I and III is interchanged in this case. It is 
also worth mentioning that in the non-interacting parts 
of the wire the Green functions are both Galilean and 
translationally invariant, depending on coordinates and 
times via £ — £' only. 



for the plasmon intensity (i — 1,2 for the left and right 
contact, respectively). The corresponding transmission 
coefficients are 71 = 1 — 7£j. It is instructive in this con- 
text to compare our present approach with that devel- 
oped in Ref. |20 . where we analyzed the tunneling density 
of states and focussed on the case of smooth variation 
of g(x) in the contact regions. As we are going to show, 
the method of Ref. ]2fJ can be generalized to the case of 
arbitrary contacts (this was briefly discussed at the end 
of Ref. I20f ) and is also useful for the analysis of the elec- 
tron distribution function. Within that approach, the 
propagator of bosons is calculated in momentum space 
(rather than in real space as in the above calculation). 
The Keldysh component of the propagator is then char- 
acterized by distribution function functions B^ 3 (to) and 
B v (uj) associated with poles at q = 7juj/v and q = -quj/u 
and describing "ghosts" (free electron-hole pairs) and 
plasmons, respectively^. While the distribution func- 
tion of ghosts is simply determined by that of incoming 
electrons, the plasmons experience in general reflection 
at the boundaries. We have for the left boundary (see 
Fig.© 



Br 



Tl x B 



L ' 



nout 



■TiBl, 
(35) 

and similarly at the right boundary. Here we have in- 
troduced the notation B™ for plasmon distributions in 
the interacting region of the wire and B° nt for out-going 
channels. Solving these equations, we find the plasmon 
distribution functions of right-movers in the interacting 
part of the wire, as well as in the outgoing channel (in 
the right lead): 



Bl = 



Br = 



71 4°) 

1-KiRa R 

T i r 2 p(0) 

i-ni%2 R 



i 



T 2^1 d(0) 

+ T 2 - 2J[T 2 



(36) 



Bf>. (37) 



The corresponding results for left movers are obtained by 
exchanging the indices R<->L and l<-»-2. 




FIG. 2: Distribution functions of plasmons B v in different 
parts of the wire. The distributions of incoming plasmons are 
determined by respective leads, B'™ — Bf, 



IV. ARBITRARY BOUNDARIES 

We turn now to generalization of these results for the 
case of an arbitrary shape of g{x) in the contact region 
between the interacting part of the wire and the non- 
interacting leads. The contact regions are in general char- 
acterized by some reflection coefficients Ti(u>) for the plas- 
mon amplitude, yielding reflection coefficients IZi — \r,i\ 2 



The method of Ref. ]20| allows us to express the expo- 
nents Tfj- in terms of these distribution functions. For 
the interacting part of the wire, we get 



Jo w 



(BX-4 0) )(l-cosu,t) 



7 (B^ + B^) (1- cos ujt)±i sin ujt 



(38) 
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The result for the tunneling into the non-interacting re- 
gion III of Fig. Q] can be obtained from Eq. ([55]) by using 
the distribution functions B° ut corresponding to this re- 
gion and replacing the interaction constant 7 by zero, 



Tr 



r^(5r-<xi-co S ^) 

Jo u 

/ — 11(1- cos <Jt)[B^\u)-B^(u)], (39) 
Jo u 



where TZ is the total reflection coefficient on a double-step 
structure, TZ = 1 - 7^/(1 - KiK 2 ). 

For the case of sharp boundaries the reflection and 
transmission coefficients are given by the Fresnel law, 
TZ h2 = (1 - K) 2 /{1 + Kf and T h2 = iK/(l + K) 2 , 
so that Eqs. (|38| and (|39f reduce to the earlier results 
(|3"!?|) . The total reflection and transmission coeffi- 
cients TZ and T take in this case the values (|34| (which 
explains the notations introduced there). Clearly, the 
general formulas (f38|) and (f39"|) can also be obtained in 
the framework of a real-space calculation that was pre- 
sented above for sharp boundaries. To do this, one has to 
modify the boundary conditions for the Green function in 
QIj in Eq. (|22|) by including the appropriate reflection and 
transmission amplitudes r,(o;) and U(u>) at two bound- 
aries and then proceeding in the same way as in course 
of the derivation of Eqs. ([26]) and ([32]) . The two methods 
(real space and k space) are thus in full agreement with 
each other. 

The formal results obtained thus far can be imple- 
mented to obtain physical observables. Consider first 
the non-interacting part of the setup, region III of Fig.Q] 
The effect of the interaction there amounts to modifi- 
cation of the distribution function of outgoing particles 
(right-movers), which has (in time domain) the form 



n R (t) = nji,o(t)e^W , 



(40) 



where Tr is given by Eq. 



This yields 



n R (t) 



I -ieVt/2 



Tr 



R 



sinh 7rTfjt + iO 



sinh tiTjJ, + iO 



(41) 



The way in which the electron distribution function is 
modified depends on the kinetics of the plasmons inside 
the interacting region. For adiabatic switching of inter- 
action, there is essentially no plasmon scattering. There- 
fore, the total reflection coefficient TZ and, consequently, 
the exponent J-r in the region III vanish. In this case 
the fermions retain their distribution function: the right- 
movers going out into the right lead have the same dis- 
tribution as the right-movers injected into the interact- 
ing region from the left lead. (The same applies to the 
left-movers, of course.) Let us now discuss the opposite 
limit of strong reflection, TZ — ► 1. For a structure with 
a sharp boundary, this is the case provided the interac- 
tion is strong, K — > 0. Alternatively, this limit may be 
realized if the boundary regions are sufficiently extended 
and characterized by random K (x) such that plasmons 
with relevant frequencies are localized. Regardless of the 
cause, in the limit TZ — ► 1 the left- and right-moving elec- 
trons exchange their distribution functions, except for 
keeping their total flux (i.e. the chemical potential). 

Next, we consider the interacting part of the wire. An- 
alyzing the result (|38|) . we see that two terms in square 
brackets have distinctly different physical origin. The 
second term, which is proportional to the local strength 
of the interaction 7 at the measurement point is respon- 
sible for creation of the zero-bias anomaly (ZBA) as well 
as for its dephasing smearing, with the non-equilibrium 
dephasing rate2£ 



7T7 



(l-ftiXl + fta) 



1 -KiK 2 
(l+KiXl-K^ 
l-TZiTZ 2 



Tr 



(42) 



On the other hand, the first term in the integrand of ( 
which is governed by the difference between the incom- 
ing and local distribution of plasmons, is fully analogous 
to the expression for J-^ in the non-interacting region, 
Eq. (f3"9"| . and describes the modification of the distribu- 
tion function inside the wire, 



n v (t) = n, h0 (t) exp j- jf° ^[B£M - ^ 0) (^)](1 - coswf) j 



1 



2?r t + iO 



exp 



dio 



— [B™{u)- l](l-coswf) 



(43) 



As is clear from Eq. (|4"3"|) , the "ghost" term with ' 
essentially serves to cancel the bare distribution function 
n v> o, so that the distribution function n(t) is determined 



only by the plasmonic distribution B™(lu) in the wire. 
This is in fact a manifestation of a general relation be- 
tween the functional and full bosonization approaches, as 
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will be discussed in detail elsewhere 2 ^. 
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FIG. 3: Total distribution functions of electrons in the left 
and right leads for the LL interaction parameters K = 1 (no 
interaction) and K = 0.2 (with sharp boundaries). Temper- 
atures of the leads are Tl = 0.2 and Tr = 0.001; the bias 
voltage is eV = 0.25. 

Fourier transformation of our results into the energy 
representation can be done numerically (for analytic cal- 
culation at equilibrium see Appendix |B| ; representative 
results are shown in Figs. [3j [4] In Fig. [3] we present dis- 
tribution functions for non-interacting parts of the wire. 
Temperatures are set to Tl — 0.2 and Tr = 0.001 (in 
arbitrary units), the applied voltage is eV = 0.25, and a 
sharp variation of the interaction at the boundaries (as 
in Sec. IIIip is assumed. The distribution function of free 
fermions (K = 1), plotted by a dashed line, is the same 
on both ends of the wire. For interacting electrons (we 
choose the interaction parameter to be K = 0.2, which is 
in the range of characteristic values reported for carbon 
nanotubes, see, e.g., Ref. Ifa) the distribution functions in 
two leads are different. In particular, the distribution 
function in the left lead (region I in Fig. [1]) has a sharp 
edge at the energy e = fj, + eV/2, which corresponds to 
cold right-moving electrons. In the right lead (region 
III), this edge is broadened due to interaction with hot 



left-moving electrons. The situation is opposite for left- 
moving particles. The distribution in the right lead has 
a broad edge at e = /i — eV/2 that corresponds to hot 
left-moving electrons. Due to interaction inside the wire 
this edge in the region I sharpens. 

In Fig. H] we present the results for the distribution 
functions of left- and right-moving quasiparticles in the 
central (interacting) part of the wire, Eq. (|43[) . For 
K = 0.2 the plasmon reflection at the boundaries is 
strong. In a symmetric structure this leads to almost 
equal distribution functions of both types of carriers in- 
side the wire. 

In the upper panel of Fig. [5] we show the results for 
TDOS for K = 0.8. The minima of TDOS are reached 
at energies e = /j,±eV/2. The broadening of the ZBA dips 
has two origins: smearing of the distribution function and 
dephasing. While the dephasing broadening [cf. second 
term in Eq. (|38|) ] is the same for both chiral branches, 
the distribution functions [cf. first term in Eq. (|38j) ] are 
in general different. A deeper minimum at e = [i + eV/2 
reflects the fact that right-moving electrons in the wire 
have a much narrower distribution function. This is be- 
cause at K = 0.8 the energy relaxation at the bound- 
aries is quite weak, so that the distribution functions of 
cold right-movers and hot left-movers are only slightly 
modified. The situation is different for K = 0.2, when 
distribution functions n/j and hl are nearly identical (up 
to a shift by eV), see Fig. [4] As a result, the structure 
of the TDOS also becomes symmetric. In fact, for the 
chosen value of the voltage, two broad ZBA dips merge 
together. 

V. THERMAL CONDUCTIVITY AND 
ELECTRONIC DISTRIBUTION FUNCTION 

We discuss now a relation between our results for the 
electron distribution function and previous findings on 
the electric and thermal conductance of a LL wire. In 
the absence of backscattering the number of left and right 
moving particles is separately preserved. As a result, the 
electric current is linear in the voltage V , 

e 2 

I = ev(N R - N L ) = —V (44) 

with unrenormalized Landauer conductance G = e 2 /h, 
Ref. [26|. In our formalism, this relation immediately fol- 
lows from Eq. (j40|) and the condition T v {t — > 0) — > 0. 
This ensures that the modification of the distribution 
function of right (or left) movers by a spatially varying 
interaction does not affect the integral of the distribution 
function over energy, i.e. the total number of carriers of 
each type. 

We turn now to the thermal conductance. The en- 
ergy current is easily found from the Green functions of 
electrons in non-interacting parts of the wire, 

I E =vd t [G<{t,t')-G<{t,t')]\ t=v , (45) 
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FIG. 4: Distribution functions of left- and right- movers, 
Eq. (143 |l . in the interacting part of the wire. All parameters 
are the same as in Fig. [3] 

which can be rewritten in terms of the electron distribu- 
tion functions, 

Ie= r ^e[n L (e)-n R (e)]. (46) 

Substituting the result (|40|) . (|39|) for the distribution 
functions, we get the expression of the thermal current in 
terms of distribution functions of incoming electron-hole 
pairs, 

i r°° 

Ie = ^J o <kvu,T(u,) [4 0) ( W ) - B£> (w)]. (47) 

According to (|47p , the thermal conductance is affected by 
the interaction [through the reflection coefficient T(u>)], 
as was first found in Ref. H3. Note that due to the 
particle-hole symmetry of LL model, the applied volt- 
age drops out of Eq. (|4"T|) . For the case of sufficiently 
sharp boundaries, when T(uj) can be considered as uj- 
independent for relevant frequencies, Eq. (|46j) reduces to 

I E = ^T{Tl-Tl). (48) 

Deviation of the transmission coefficient T(uS) from unity 
leads to the violation of the Wiedemann- Franz law 27 . As 




Q| 1 1 1 1 1 1 1 
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8 

FIG. 5: TDOS u(e) (normalized to its non-interacting value 
2^o) in the interacting region for K = 0.8, 0.4, and 0.2. 
The temperatures of leads and the voltage are the same as 
in Fig. H 



is seen from our analysis, this deviation is a manifesta- 
tion of a microscopic phenomenon: energy relaxation of 
electrons due to non-uniform interaction. 

The heat current ([47)) can be equivalently represented 
in terms of plasmonic distributions in the wire 

I E = ^J°°dwu J [Bl(u } )-Bl(u;)]. (49) 

This implies that the presentation of the heat current 
in the form (|46|) is also valid in the interacting part of 
the wire, with the electronic distribution functions n n (e) 
given by (|43|) . Thus, also in the interacting part of the 
wire, the energy current can be understood as carried by 
properly defined quasi-particle excitation. This is a re- 
markable result, which demonstrates that the concept of 
fermionic quasiparticles remains meaningful in a strongly 
interacting ID system (LL) despite its non-Fermi-liquid 
features. 
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VI. SUMMARY AND OUTLOOK 

To summarize, we have developed a theory of tunnel- 
ing spectroscopy of LL conductor connected to reservoirs 
away from equilibrium. In the specific setup considered 
here, each branch originates from a source which is at 
equilibrium. However, the right and the left sources have 
different temperatures and different chemical potentials. 
We have modeled the system as a LL with spatially non- 
uniform interaction, and calculated the single-electron 
Green functions that carry information about the 
TDOS and the fermionic distribution functions in differ- 
ent parts of the wire. The interaction affects the tun- 
neling characteristics in three distinct ways. First, it in- 
duces a power-law ZBA in the TDOS v(e) (with two dips 
split by the voltage) in the interacting part of the wire. 
Second, it leads to broadening of ZBA singularities due 
to dephasing, with the dephasing rate governed by the 
interaction strength and the plasmon distribution inside 
the wire. Both the ZBA and the dephasing effects are 
encoded in the second term of Eq. . 

The third effect of the interaction — which is specif- 
ically at the focus of the present work — is the inelas- 
tic scattering of electrons, leading to their redistribution 
over energies. This effect takes place in those regions 
where the interaction strength varies in space (near the 
wire boundaries in our model), inducing backscattering 
of plasmons (but not of electrons). This leads to relax- 
ation of the electron distribution functions: left and right 
moving fermions "partly exchange" their distributions, 
see Eqs. (@3| and Figs. 0] [5] For slowly varying 

interaction, when the plasmons with relevant frequencies 
go through essentially without reflection, the energy re- 
laxation of electrons is negligible. In the opposite limit, 
when the plasmons are almost entirely reflected (due to 
strong and sharply switched interaction or, else, due to 
disordered boundary regions inducing the plasmon local- 
ization), the left- and right-movers essentially exchange 
their distribution functions (but not their total density). 
We have also discussed a connection between these re- 
sults and earlier findings on the thermal conductivity of 
LL structures. 

Our results are important for the analysis of TS exper- 
iments on strongly correlated ID structures (in particu- 
lar, carbon nanotubes?) out of equilibrium. In this con- 
nection, let us emphasize the following important point. 
What can actually be measured in experiment are Green 
functions, G> and G< . The TDOS v(e) in the interact- 
ing part of the wire, as well as the distribution function 
n(e) in the non-interacting regions are related to G > and 
G < in a simple way. On the other hand, in order to ex- 
tract the distributions n/j(e) and n_t(e) from G^ in the 
interacting part of the wire, a non-trivial deconvolution 
procedure is necessary. The broadening of (split) Fermi- 
edge structures in G< in the interacting part of the wire 
is governed by both the distribution function and the de- 
phasing. The dephasing contributes to the smearing of 
Fermi-edge singularities also in higher-dimensional (diffu- 



sive) systems 2 ^, and should be taken into account for the 
accurate interpretation of corresponding experiments^ 2 ^. 
In the ID case the role of dephasing becomes particularly 
dramatic (if the interaction is sufficiently strong). This 
is very well illustrated by Fig. [5] two Fermi-edge singu- 
larities almost (middle panel) or even completely (lower 
panel) merge, despite the fact that the Fermi edges in 
the distribution functions remain well separated (Fig. [5]). 

A comment of a more general nature is in order here. 
Our results illustrate the fact that there is no unique an- 
swer to the question: "How much is a LL different from 
a Fermi liquid?" On one hand, the strong, power-law 
ZBA in TDOS of a LL clearly distinguishes it from the 
Fermi liquid. In more formal terms, the single-particle 
residue Z, which is finite in the Fermi liquid, vanishes in 
a power-law fashion at the Fermi level of the LL. Also 
the dephasing rate determining the broadening of ZBA, 
Eq. (|42[) . is linear in temperature, contrary to the Fermi- 
liquid T 2 behavior. One could think that it makes little 
sense to speak about fermionic excitations in this situa- 
tion, but this is not the case. First, the power-law van- 
ishing of TDOS has little importance (like the value of Z 
in the Fermi liquid) for kinetic properties of the system. 
Second, the dephasing rate (|42[) is governed by processes 
with zero energy transfer and do not lead to any en- 
ergy relaxation. As a result, the distribution function of 
fermionic excitations, n^ie), is a fully meaningful concept 
even in the case of a strong interaction. It stays preserved 
as long as the interaction is spatially constant (or varies 
adiabatically slow with x). Furthermore, both the charge 
and the energy current in the interacting part of the wire 
can be understood as carried by these fermionic quasi- 
particles. From this point of view, the LL is a perfect 
Fermi liquid. 

We conclude the paper by reviewing some future re- 
search prospects; the work in those directions is currently 
underway. First, one may consider a more general non- 
equilibrium situation where the distribution functions 
"injected" into the interacting part of the wire are of non- 
equilibrium (e.g., double-step) form by themselve o 20 ' 24 , 
see setups b, c in Fig. I of Ref. 20; the first of these 
setups is close to the experimental situation of Ref. H. 
This requires a generalization of the bosonization tech- 
nique that will be presented elsewhere 2 ^. Second, it is 
interesting to study correlations between outgoing left- 
and right-movers. In a general situation, one finds that 
their density matrices are not decoupled, i.e. they are 
entangled, which manifests itself, in particular, in cur- 
rent cross-correlations. Third, one may study the effect 
of a random variation of the interaction strength K(x) in 
the wire. If the wire is sufficiently long, plasmons with 
not too low frequencies get localized. Using our general 
results, one concludes that in the left (right) half of the 
wire both distributions hr, are determined by that 
of the left (respectively, right) reservoir, with a transi- 
tion region which extends over the localization length of 
the middle section. To refine this picture, one has to 
include into consideration also plasmons with low fre- 
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quencies (that remain delocalized). Also, including the 
spectral curvature will induce plasmon decay processes. 
(In the context of thermal conductivity, this physics was 
discussed in Ref. [2?].) Finally, our results can be general- 
ized to the case of chiral LL, where both branches move in 
the same direction, which is the situation characteristic 
for quantum-Hall edge-state devices^. 
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where the subscripts "tp" and "w" refer to the tun- 
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tation. If electron tunneling is local in space, we have 
T(y,y') — T6(y — y')5(y — x), where x is a position of 
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the probe, 
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APPENDIX A: MEASUREMENT OF GREEN 
FUNCTIONS G% 



Gf p {e,x,x) = 27iwtp(e)n t p(e), 

G>{e,x,x) = -27rii* p (e)[l-nt p (e)]. (A2) 



The tunneling current between a probe and a quantum 
wire can be expressed in terms of the functions G^ as 

I(U) = J dydy'\T y ^\ 2 J * 

x G<(e-eU,y,y')G>(e,y',y) 



GfM- eU,y,y')G<(e,y',y) 



(Al) 



Differentiating the tunneling current with respect to 
voltage and substituting Eq. (|A2[) into Eq. (|A1[) . one finds 



dl_ 

dU 



-2i\T\- 



de 



-27iw tp (e - eC/)i/ w (e) 



dutpje - eU) 
de 

dn tp (e - eU) 
de 



n tp (e - e£/)G> (e, x, x) + (1 - n tp (e - eU))G< (e, x, x) 



(A3) 



For a LL wire the Green functions G w and the TDOS v w 
represent a sum of contributions of both chiral branches. 
If the density of states in the tunneling probe (ftt>) IS a 
constant (as in a normal metal) , the first term in Eq. (|A3|) 
drops out. In this case the result is proportional to the 
TDOS in the wire. Assuming that the tunneling probe 
is at zero temperature, one then finds 

^-=4n\T\^ w (eU). (A4) 
On the other hand, if the density of states in the tun- 



neling probe is strongly energy dependent (as for super- 
conducting electrodes), the first term in Eq. (|A3|) sur- 
vives. Unlike TDOS (which is determined by the differ- 
ence G^ — G^), this term contains also the information 
about G^+G^. Therefore, measurement of the tunneling 
current with two different types of tunneling probes (nor- 
mal and superconducting) allows one to find functions 
G> and G^ separately. The idea to use superconducting 
electrodes for the tunneling spectroscopy was introduced 
m Ref.0 and more recently employed in Ref. [8|. 
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APPENDIX B: GREEN FUNCTIONS G< AT 
THERMAL EQUILIBRIUM 



This yields the following asymptotic behavior of the 
Green function at low temperatures (|e| 3> T), 



At thermal equilibrium the Green functions in the en- 
ergy domain can be calculated explicitly. Using Eq. 
and Bl = Bl = B% =B% = coth we find the expo- 



nent (t) for the Green functions in interacting part of 
the wire, 



iAsmhnT(t - i/A) 



(Bl) 



where we drop the chirality index 77, as it is immaterial for 
x = x' in equilibrium. Using Eq. (|18[) and performing a 
Fourier transform from the time into the energy domain, 
one finds 



G' 



(ttT) 



1+7 



2-Kv{iA)~< J_ 



dte 1 



sinh i+7 nT(t-i/A) 



After calculating an auxiliary integral 

P izt 

dt 



(B2) 



sinh 1+7 (i-i0) 



^1+727 
' — ; r< 

r(i + 7) 



= 7T2/2 



r[(l + 7 + zz)/2] 



one obtains 



G > (e) = -- 



2"' 



ttT 



r[(l +7 + zz)/2] , 
(B3) 



27rvr(i + 7) V a 

where z = e/wT. Similarly, one finds the function G < , 



27r^r(i + 7) \ a J 



r[(l + 7 + w)/2] 

(B4) 



G>(e) 



1 *(z-\z\)/2 ( l£i 

A 



«r(i + 7) 



(B5) 



and high temperatures (|e| <C T), 



G>(e) 



Using Eqs. (J3T|), (jB3|, and (pi]) , one obtains TDOS at 
equilibrium, 



^(e,T) 



2T- 1 /ttT 



7r 2 t>r(l + 7) V A 
r[(l + 7 + iz)/2] 



7TZ 

cosh-^-. (B7) 



Equation (jB7j) describes the well-known ZBA in TDOS, 
u(e) cx |e| 7 , smeared at the scale e ~ 2itT(1 + 7). This 
smearing results from a combined effect of (i) the ther- 
mal broadening of the distribution function and (ii) the 
dephasing rate2£ l/r^ = 2njT. 

It is straightforward to check that the Fermi-Dirac dis- 
tribution function is recovered from the ratio 



in agreement with the fluctuation-dissipation theorem. 
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